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METHOD FOR EVALUATING ELECTRICAL FffiLD LEVELS 
IN THE NEAR-ZONE OF TRANSMITTING ANTENNAS 

5 Description 

This invention relates to tools for planning radio carrier telecommunications 
systems and particularly relates to a method for evaluating electrical field levels 
in near-zones of transmitting antennas. 

10 

The increase of radio base stations (RBSs) for mobile communications in highly 
populated areas has caused grovdng anxiety in people living close to said 
stations^ worried about biological effects from electromagnetic pollution on 
the human body. As a consequence^ governments in various countries have 
15 established stringent exposure limits for the level of electromagnetic fields 
radiated in the radio-frequency range, in zones where the population may be 
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exposed for long periods of time. Mobile telephony operators must produce 
more detailed documentation on the compliance of their radio base stations 
to prescribed limits. Consequently, the use of more accurate and reliable 
planning tools is necessary to avoid costly and difBcult interventions on devices 
5 and systems after construction and installation. 

The electrical field radiated by transmitting antennas, particularly radio base 
station linear arrange antennas, is currendy calculated under a far field 
hypothesis, expressing an antenna gain in spherical co-ordinates G(6,(|>) as the 
10 product of radiation patterns in the horizontal and vertical planes. 

The far field hypothesis is verified when the distance from the antenna exceeds 
— ^, where is the maximum extension of the antenna and X is the wave- 

length. Considering base radio station antennas and cellular telephony system 
15 frequency bands, this distance can correspond to more than 40 metres, while 
in urban areas base radio stations are usually positioned at less than 10 metres 
from inhabited buildings. Consequentiy, the use of more accurate methods for 
evaluating electrical field levels in the near-zone of radio-frequency sources, 
particularly linear array antennas, vridely used in cellular telephony systems, 
20 is necessary. 

A large number of numerical methods and software tools for the evaluation of 
the electromagnetic field radiated by antennas are currentiy available. For 
example, the well-known software NEC-2, described in 'TMumerical Electromag- 
25 netic Code - Method of Moment", by Burk G.J. and Poggio A.J., Lawrence 
Livermore National Laboratory, January 1981, calculates the electromagnetic 
field very carefully, both in the far-field and in the near-zone of the transmitting 
antenna, but requires a good knowledge of wire-grid technique for input file 
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implementation^ a precise mechanical quotation of all parts of the antenna and 
extensive computation time. In practice, however the knowledge of the 
antenna is limited to its physical dimensions, its gain and its radiation patterns 
in the vertical and horizontal planes. 

The method for evaluating electrical field levels in the near-zone of transmitting 
antennas which is the object of the invention overcomes said shortcomings and 
solves the described technical problems by providing an accurate evaluation 
while requiring a limited number of input data and reduced processing time 
by data processing units used for calculating. 

Spedficall}^ the object of the invention is a method for evaluating the electrical 
field levels in the near-zone of transmitting antennas, as described in the 
characterising part of claim 1. 

Additional characteristics and advantages of the invention will now be de- 
scribed, by way of example only, with reference to the accompanying drawings 
wherein: 

Fig. 1 schematically shows a linear array antenna arranged along the 
z axis of a Cartesian reference system in symmetric position with 
respect to axis x; 

Fig, 2 represents a thin wire of arbitrary length L arranged on the z 
axis of a Cartesian reference systems in symmetric position with 
respect to axis x. 

As previously mentioned, a simple method for estimating the electrical field 
radiated by a radio base linear array antenna is based on the hypothesis that 
the antenna gain G(0,(})) can be calculated as the product of radiation patterns 
in the vertical and horizontal planes. These patterns can be determined by 
means of measurements made in the first phases of the procedure, along with 
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the physical dimensional measturements on the antenna. Alternatively, the 
patterns provided by the manufacturer can be used. 

According to the method of the invention, the antenna gain G(6,4>) is calculated 
in a slightly different way by introducing the angle t|r, i.e. : 



G(e,4)) = GMDv(e)D„(!|r) 



(1) 



where: 



is the maximum gain of the antenna 
Dv(6) is the radiation pattern in the vertical plane (<() = 0) 
D„(i|;) is the radiation pattern on a conical surface generated by the 
rotation on axis z in the maximum radiation direction 
angle function of 6, ^ and x, results from: 



*(e,<l>,t) = ct)sin(P(e,T)), where p(e,T)- 



ic+2t 

7t 



7C-2T 



7c/2+t<eiic 



(2) 



where x is the down-tilt of the radiation beam with respect to the horizontal 
plane. 



The angle i|f is introduced to properly match the secondary lobes of the radia- 
tion pattern. Its analytical expression derives from setting the following 
constraints: 

The direction (6,(j)) does not depend on the angle (|) for 6 = 0, Ji. This 
is because per 0 = 0 the direction (0,<|)) is always the direction +2, 
while for 0 = n the direction (0,<|>) is always direction -z. Conse- 
quently, when 0 = 0,71 the value of G(0,<t)) must not depend on the 
angle 4> and, specifically, must be equal to the value of the vertical 
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pattern on 0 - 0,n: because Dh(0) = 1, the result is obtained by 
setting !|i = 0; 

for 6 s= ic/2 + T, the function 6(8,4)) must coincide with the radiation 
pattern in the horizontal plane. The result Is obtained by setting i|r 
= (j) because in 6 = n/2 + t the function D^CG) has unitary value. 

As known, the electrical field £(r,d,(j)) in the far-zone is given by: 



where is the radiated powe^ r is the distance from the electrical centre of 
the source and p is the polarisation vector. 

We will now suppose that the separation of the effects in the vertical plane and 
in the horizontal plane, as expressed by equation (1), is still valid in the near- 

21^ 

zone field, particularly in the zone 2 to 3 A and , where is the maximum 

extension of the antenna. Because the physical width of a linear array antenna 
is much smaller than its height, the horizontal pattern contribution D„(i|r) to 
the electrical field in the near-zone iield is approximately equal to that in a far- 
zone field. 

According to this hypothesis, the near-zone electrical field can be expressed 
as: 




(3) 



r 



|E(r,e,(|))| « v^PRD„(t|F)|E(r,e)| 



(4) 



where £(r,6) is an appropriate vectorial function which is calculated as de- 



wo 02/095428 



PCT/EP02/056S7 



6 

scribed below. 

Equation 4 can consequently be used to calculate the electrical field in a point 
P - (r^6,(t)) of the space as a product of two contributions: the first of which 
5 is due to the electrical field in the horizontal plane of the antenna D||(i|r) and 
the second of which is due to the electrical field in the vertical plane of tiie 
antenna Dv(6). 

The second contribution is obtained by calculating the electrical field in 
10 projection Pj = (r,e,0) of point P = Cr,e,<l)) on the vertical plane x, z of the 
antenna, as shown in Fig. 1. In this figure, the linear array antenna, schemati- 
cally iUustrated, is arranged along the axis z of a Cartesian reference system 
X, y, z in a symmetric position with respect to axis x. Point ? identified by 
spherical co-ordinates r,6,(|), is the point in space where the electrical field 
15 intensity is to be evaluated. Point P, is the projection of point P on the vertical 
plane of the antenna (plane x, z) and is obtained by turning point P, around 
axis z, by an angle which is equal to the value of the spherical co-ordinate 4>, 

To calculate the electrical field in Pp we introduce the concept of "Equivalent 
20 Current Distribution", abbreviated as ECD, i.e. the distribution of current vdiidi 
generates the same radiative effects of a certain source in the near-zone field 
and in the far-zone field. 

In the case of a linear array antenna, which is verticaUy polarised, the ECD can 
25 be conveniently defined as a one-dimensional cwrent along a thin wire of 
arbitrary length L, as shown along axis z in symmetrical position with respect 
to axis x in the Cartesian reference system of Fig. 2. 

In the figure, I(z) is the electrical current running in the wire, which is sup- 
30 posed polarised in direction z (i.e. I(z) = I(z) i). Vector t identifies a generic 
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current element on the thin wire. Co-ordinates 0^ ()>4 = 0 are the spherical 
co-ordinates of point Px in a Cartesian reference system whose origin corre- 
sponds to the current element identified by vector j^. 

ECD current is obtained from the knowledge of the radiation pattern in the 
vertical plane of the antenna: from the radiation theory, we can write the 
radiation pattern in terms of voltage V(e) in the vertical plane z as: 

V(e) = yDy(0) = Csine f I(2)exp027C2cos8)d2 (5) 

-1/2 

where C is a multiplying constant and I(z) is the current generated in the wire. 

By defining: 
u = (cosOyX, 

Viu) = V(u)/^l-au)^ 

and resolving (5) with respect to I(z), i.e. to ECD, we obtain: 

♦1/2 

1(2) = I V(u)exp(-j27czu)du = S{^(u)> (6) 

i.e., ECD is the Fourier transform of the previously defined radiation pattern 

V (u) . It can be calculated with a standard Fourier Fast lY'ansform (FFT), by 

using the Nyquist sampling theorem and estimating the physical extension of 
the ECD as x times the physical height of the antenna. The resulting ECD is 
consequently discrete, i.e. formed of a certain number of current samples I„. 
Nimieric results demonstrate that % values in the range from 2 to 3 are accept- 
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able for linear arrays. 

Alternatively^ the ECD can be obtained with the Woodward-Lawson sampling 
method (for example, see article by G. A. Somers "A proof of the Woodward- 
Lawson sampling method for a finite linear array**. Radio Science, VoL 28, No. 
4, pp. 481-485, July-August 1993). 

Hie ECD determined in this way is a discrete current distribution, therefore the 
vectorial function £(r,6), required to compute the electrical field in the vertical 
plane according to (4), is given by: 



where is the characteristic impedance of vacuum, equal to 377 Q, X is the 
wavelength, ko ^ 2n/X is the constant propagation in vacuum, I„ is the nth ECD 
sample, and d and 6^, as previously mentioned, are the spherical co-ordinates 
of point P^in which to compute the electrical field in the Cartesian reference 
centred in the nth current sample as shown in Fig. 2. 

Naturally, numerous changes can be implemented to the construction and 
embodiments of the invention herein envisaged without departing from the 
scope of tiie present invention, as defined by the following claims. 



expC-jnkpd) 
d 




(7) 
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daims 



1. Method for evaluating electrical field levels | E (r, e, (J)) | in a point P(i;6,<t)) 
in the near-zone of transmitting antennas, particularly linear array anten- 
nas located in the origin of a Cartesian reference system identified by the 
orthogonal co-ordinates (x, y, z) and the spherical co-ordinates te6,4>)« 
in which: 

the power (Pr) of the radio-firequency transmitted by the antenna is 
established; 

the radiation pattern (Dy) of the antenna in the vertical plane is 
determined (<|> » O); 

the radiation pattern (Dh) of the antenna in the horizontal plane is 
determined; 

the physical height of the antenna is measured; 
characterised in that said signal power (Pg) and said radiation patterns 
(Dvj D„) are used as inputs for data processing means to determine the 
electrical field level according to the following expression: 



where: D„(i|/) is the radiation pattern on a conical surface generated by 
the rotation about axis z in the maximtmi radiation direction; angle iff, 
function of 6, 4) and is given by: 



|E(r,e,<|))| = ^PgD„(i|r)|E(r,0)| 



6 



^M,x) = 4)sin(p(e,T)), wherep(0,T) = 



ic+2t 



(ic-e) 



1C-2T 
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T being the down-tilt of the radiation beam with respect to the horizontal 
plane (x, y); and £(r,6) is a vectorial function resulting from: 

where Zq is the characteristic impedance of vacutmi, A is the wavelength, 
ko is the propagation constant in vacuimi, I„ is the nth sample of one- 
dimension current (ECD) along a thin wire of arbitrary length (L), which 
generates the same radiative effects of the antenna both in the near-zone 
and the far-zone, and d and 6^ are the spherical co-ordinates of the 
projection = (r,e,0) of said point P(r,e,<I>) in the vertical plane (3C, z). 

2. Method according to claim 1, characterised in that the following Fourier 
transform 9(9^ (u)} is used to determine said one-dimensional current 
(ECD) by means of data processing means: 

I(z) = j V'MexpC-j^TczuDdu = 9{^(u)} 
where u = -1/2 

cos(e)/X and ^(u) = V(u)//l -(Au)^, V(u) being the voltage radiation 

pattern in the vertical plane (x, z), resulting from radiation theory of said 
radiation pattern (Dy). 

3. Method according to claim 2, characterised in that said Fourier transform 
is calculated by said data processing means with a standard Fourier I^t 
TY-ansform (FFT), by using the Nyquist sampling theorem and estimating 
the physical extension of said ohe-dimensiohal ciurent (ECD) as % times 
said physical height of the antenna, being x compressed in the range from 
2 to 3. 
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4. Method according to claim 2, characterised in that said one-dunensional 
current (ECD) is calculated by data processing means according to 
Woodward-Lawson sampling method. 
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P = (r,e,0) 

p, = (r,ao) 




Fig. 1 



wo 02/095428 



2/2 



PCT/EP02/05657 




0d = 0 = 



-L/2 Y 



Fig. 2 



IRNATIONAL SEARCH REPORT 



itlonaJ Application No 

PCT/EP 02/05657 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 7 G01R29/10 



AccordinQ to International Patent Classllicatlon (IPC) or to both natlortal classification and IPC 
B. FIELDS SEARCHED 



Minimum documentation searched (dassillcallon syslem toilowed by ciassiticalion symbols) 

IPC 7 601R 



Oocumentailon searched other than minimum documentation lo the extent thai such documents are Included In the fields searched 



Electronic data base consulted during the international search (name ol data base 

EPO-Internal 



and. where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category * Citation of document, with indication, where approprfaie, of the relevanl passages 



Retevant to clabn Na 



MCLEAN J ET AL: "Measurement of the 
antenna factor of diverse antennas using a 
modified standard site method" 
ELECTROMAGNETIC COMPATIBILITY, 1999 IEEE 
INTERNATIONAL SYMPOSIUM ON SEATTLE, WA, 
USA 2-6 AUG, 1999, PISCATAWAY, NJ, 
USA, IEEE, US. 

2 August 1999 (1999-08-02), pages 
595-600. XP010360732 
ISBN: 0-7803-5057-X 

page 595, left-hand column, paragraph 1 
-page 597, left-hand column, paragraph 1; 
figure 1 



I I Further documents ere listed In the continuation of box a 



I I Patent family members are listed bi 



* Special categories of died documents : 

•A* document defining the general stale ol the art which Is not 

considered to be of panicular relevance 
'E' earlier document but published on or after the Intemallonal 

filing date 

•L' document which may throw doubts on priority clalm(s) or 
which is cited to establish the pubficailon date of another 
citation or other special reason (as specified) 

'O* document referring to an oral disdosura. use, exhibition or 
other means 

*P* document published prior to the international fiflng date but 
later than the priority date claimed 



T later document published aflef Ihe bitematkmal nDng dale 
or priority dale and not tn conflid wRh the application but 
died U) understand the prfndple or theory underlying the 
invention 

*X' document of panicular relevance: the daimad Invention 
cannot be considered novel or cannot be considered lo 
involve an inventive step when Ihe document Is taken i 

•y document of particular relevance: the claimed invention 
cannot be considered to invoke an inventive step when the 
document Is combined with one or more other such docu- 
ments, such combination being obvious to a person sMIIed 
ki the art 

*A' document member of the same patent family 



Date of the actual completion of the Intemational search 



25 September 2002 



Date of mailing of the Intemational search report 

08/10/2002 



and mailing address of the ISA 

European Patent Office. P.B. 5616 Patenllaan 2 
NL>2280HVRisvtrQk 
TeL (♦31-70) 340-2040. Tx. 31 651 epo nl. 
Pax: (431^70) 340-3016 



Authorized officer 



Heinslus, R 



PCTASA/210 (awond ahMi) (July tflOS) 



This Page is inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 



BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the 
original documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

^ BLACK BORDERS 

IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 
^ FADED TEXT OR DRAWING 

□ BLURED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

^ COLORED OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REPERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: ^ 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning documents will not correct images 
problems checked, please do not report the 
problems to the IFW Image Problem Mailbox 



